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Abstract 

Fractional mechanics describes both conservative and non-conservative 
systems. The fractional variational principles gained importance in study- 
ing the fractional mechanics and several versions are proposed. In classical 
mechanics the equivalent Lagrangians play an important role because they 
admit the same Euler-Lagrange equations. By adding a total time deriva- 
tive of a suitable function to a given classical Lagrangian or by multiplying 
with a constant the Lagrangian we obtain the same equations of motion. In 
this study, the fractional discrete Lagrangians which differs by a fractional 
derivative are analyzed within Riemann-Liouville fractional derivatives. As 
a consequence of applying this procedure the classical results are reobtained 
as a special case. 

The fractional generalization of Faa di Bruno formula is used in order 
to obtain the concrete expression of the fractional Lagrangians which dif- 
fers from a given fractional Lagrangian by adding a fractional derivative. 
The fractional Euler-Lagrange and Hamilton equations corresponding to 
the obtained fractional Lagrangians are investigated and two examples are 
analyzed in details. 

Keywords: fractional Lagrangians, fractional calculus, fractional Riemann-Liouville 
derivative, fractional Euler-Lagrange equations, Faa di Bruno formula. 

1 Introduction 

Fractional calculus is a generalization of ordinary differentiation and integration 
to arbitrary order. The fractional derivatives are the infinitesimal generators of 
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a class of translation invariant convolution semigroups which appear universally 
as attractors. Various applications of fractional calculus are based on replacing 
the time derivative in an evolution equation with a derivative of fractional order. 
The results of several recent researchers confirm that fractional derivatives seem 
to arise for important mathematical reasons. During the last decades the frac- 
tional calculus [1-3] started to be used in various fields, e.g. engineering, physics, 
biology and many important results were reported [4-19]. A special attention 
has been devoted during the last years to the fractional variational principles 
[20-29] and their applications [30-35]. For the constrained systems [36] the frac- 
tional Lagrangian and Hamiltonian formalism are still at the beginning of their 
development. The formulation of the fractional variational principles still needs 
to be more elaborated in the future and it will have an important impact on 
the elaboration of a consistent fractional quantization method for both discrete 
and continuous systems. The fractional variational principles is deeply related to 
the fractional quantization procedure. For a given Lagrangian there are several 
proposed methods to obtain the fractional Euler-Lagrange equations and the cor- 
responding Hamiltonians. However this issue is not yet complectly clarified and 
it requires more further detailed analysis. One of the main obstacles is related to 
the fractional Leibniz rule, the fractional chain rule [37] as well as the fractional 
Taylor series [38, 39]. The second problem is related to the non- locality of the frac- 
tional Lagrangian. Recently, the Hamiltonian formalism for nonlocal Lagrangians 
was investigated in [40] and a relation of non-local theories and Ostrogradski's 
formalism [41, 42] was reported [43], [44]. 

As it is well known the classical equivalent Lagrangians represent a powerful 
tool for studying the Hamilton- Jacobi equation in C aratheodory' s formulation 
[45] as well as in generalized mechanics [46, 47]. 

For these reasons, in this study we generalized the notion of equivalent La- 
grangian for the fractional case. The fractional Euler-Lagrange equations of frac- 
tional equivalent Lagrangians are calculated and the fractional Hamiltonians are 
constructed. 

The plan of this paper is as follows: 

In Section 2, some basic formulas of the fractional calculus are briefly reviewed. 
Section 3 briefly review the fractional Hamiltonian approach of discrete systems. 
In Section 4 the fractional Lagrangians are investigated and two examples are 
investigated. Section 5 is dedicated to our conclusions. 

2 Brief overview of fractional calculus 

In this section, we formulate the problem in terms of the left and the right 
Riemann-Liouville (RL) fractional derivatives, which are defined as follows: the left Rie- 
mann -Liouville fractional derivative 
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D "W = r5^(l)"/(- T + t ^ 1 ^*- (1) 



and the right Riemann -Liouville fractional derivative 

' D ?/w-f(^)(-l)7 (T - f) ""°" l/wrfT ' 



(2) 



where the order a fulfills n — 1 < a < n and V represents the Euler's Gamma 
function. It is observed that if a becomes an integer, we recovered the usual 
definitions, namely, 

«D°/(t) = [^j fit), f(t), a = 1,2,.... (3) 

Fractional RL derivatives have many interesting properties. By direct calculation 
we observe that the RL derivative of a constant is not zero, namely 



, (t-a)- a 

r(i-a) 



a D°C = C v TVl '— (4) 



RL derivative of a power of t has the following form 

^ - r{p _ a + 1 y (5) 

for a > — 1, P > 0. Composite of fractional derivatives is given by the following 
formula 



(t-a) 



j=i l \ i a j j 



-ot-J 



(6) 



where 0<k — l<q<k,p>0 and k is a whole number. As it can be seen from 
(6) the composition is not commutative. Finally, the fractional product rule is 
given below 

jw 9 )=£ f:)(^/)(§). Pi 



3=0 



By inspection we observe that the fractional product contains infinitely many 
terms and this product is taking into account the memory. 



2.1 Fractional derivative of a composite function 

In order to find the fractional generalization of a composition function we have 
obtain first of all the most general classical counterpart. 
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Let us take an analytic function <f>(t) and f(t) = H(t — a), where H(t) is 
the Heaviside function. Using the Leibniz rule and the formula for the fractional 
differentiation of the Heaviside function we obtain 

o D?0(f) = E ( fc ) <!> ik) (t) a T>r k H(t - a) (8) 

or 

under the assumption t > a. Let us suppose that <p(t) is a composite function 

<Kt)=F(h(t)). (10) 

The k-th order derivative of 4>{t) is evaluated with the help of the Faa di 
Bruno formula [37] 

Jk k k -i fh( r )(i\\ ar 

^F(h(t)) = k\^F^\h(t))Y,J[^[^) , (ID 

where the sum J2 extends over all combinations of non-negative integer values of 
a\ , • • • , cik such that 

k 

J2ra r = k (12) 

r=l 

and 

k 

a r — m. (13) 

r=l 

Therefore, the fractional derivative of a composition function is given by 



(14) 



Mi) k '^i^ M ^m 

where the sum an d coefficients a r have the meaning explained above [2]. 

3 Fractional Lagrangian and Hamiltonian 
analysis of discrete systems 

For a given fractional Lagrangian given by 

Lftfi^aDWitUBfcfit)), p= N, (15) 
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the Euler-Lagrange equations are given as follows (see for example Refs. 
[22], [28] and the references therein) 



fe + «d^ + ^dMir 0,0<o<1 - (16) 



For simplicity,in the following we consider the following form of the fractional 
Euler-Lagrange equations 



^+' D ? a iv»^ =0 ' 0<a<l,p=l,...,JV. (17) 



dqP(t) b d a V«qP{t) 

In the following by using (17) we define the generalized momenta as (see 
Ref.[27] for more details) 

p - = a^)' p=1 '" , ' iv - (18) 

As a consequence of (17) and (18) a Hamiltonian function is defined as 

H = p apa Btq p (t)-L' f . (19) 
The canonical equations corresponding to (19) are given below 



«- = -*• a^ = " D ' 9 '' a ? = ^.o<«<i,P = i,-,w. (20) 

We mention that other interesting formulations of fractional Lagrangian and 
Hamiltonian dynamics can be found in [20, 21] and [23]. 



4 Fractional derivatives and fractional Lagrangian 
mechanics 

Let L(q p (t) , q p (t)) be a classical Lagrangian function with p = 1,---,N. Let 
L'(qP,q p (t)) = L(qP(t),qP(t)) + where F(q m (t)) is a differentiable func- 

tion and q m (t) is one of the coordinates. It is very well known that two equivalent 
classical Lagrangians admit the same Euler-Lagrange equations. In the following 
we are considering the fractional generalization of this classical results, namely 
we replace the normal derivatives by the fractional ones. 
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4.1 A particular case 

One of the possible generalizations of the classical Lagrangian L(q p (t),q p (t)) is 
given by 

L f tf(t),J%<f(t)). (21) 
If we add to (21) a term of the form t D"g(i) we obtain 

L' f = L f (qO(t), a T>W(t)) + C a B?q m (t), (22) 

where m G {1, • • • , N} is a chosen coordinate and C is a non-zero real constant. 
The fractional Euler-Lagrange equations of (22) are given by 

a ^ + .D^t^ = 0, (23) 



dcftt) 1 w<r(t) 

9L < + 0*^ + 0^ = 0, (24, 



dq™(t) b d a T>?q™(t) T(l - a) 
dL f ^ 9L f 



^) +tD? WW-°' (25) 

where a — 1 , • • • , m — 1 and 5 = m + 1 , • • • , N. 

The last term of (24) arises after taking into account (16) and (4) with C — 1. 
4.1.1 An example 

Let us consider the fractional generalization of a free Lagrangian of one degree of 
freedom L — ^- at which we add a term of the form x. At the classical level the 
Lagrangians L = ^ and L' = ^ + % are equivalent. 

One of the possible generalizations for the fractional case is given below 

L f = iaB? f ))2 +C a B?q(t), (26) 

where C is a non-zero real constant. 

By using (26) the fractional equation of motion is given bellow 

t Df (^(0) + C f0^ = °- (27) 
From (26) the expression of the fractional canonical momentum is 

Pa = a T>«q{t) + C. (28) 
The fractional canonical Hamiltonian is given by 

H f = Paa BMt)-L f = iPa ~ C)2 . (29) 
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Applying the method developed in [27] we obtain 



an f)fj 

- = a U ?q(t) ,- = ,U^, (30) 

or 

a -D?q(t) = Pa -C, tDfp a = 0. (31) 

Therefore, from (31) and (27) and taking into account (28) we conclude that 
fractional Euler-Lagrange and fractional Hamilton equations give the same evo- 
lution equation for the coordinate q(t). Besides, when the constant C becomes 
zero and a — > 1 the classical results are reobtained. 



4.2 The general case 

The next step is to add to the fractional Lagrangian Lf(q p (t),aD"q p (t)), p = 
1, • • • , N a most general term which under some limits becomes the classical total 
derivative. We analyze the fractional equivalent Lagrangians as it is given below 

L' f = L f (q<>(t), a T><?q»(t)) + a B?F(q m (t)), p = 1, • • • , n. (32) 

By using the fractional generalization of Faa di Bruno formula given in (14) we 
obtain the explicit form of (32) as follows 



L) = L f (q?(t), a BW(t)) + (t))+ 

(t-a) k - a * s) ^ u * 1 f(q m Y r) (t)\ ar 



r! 



.(33) 



The presence of the infinite higher order derivatives in (33) rises an interesting 
question regarding the form of the corresponding Euler-Lagrange equations as well 
as the form of the corresponding Hamiltonian construction. The theory described 
by (33) is non-local and involve infinity derivative of q(t). The key point is to use 
the formula (9) which allowed to represent the both parts of the Lagrangian in a 
common way. 

The next step is to write the Euler-Lagrange equations and the corresponding 
Hamiltonian. 

In this line of taught we consider that the dynamical variable q(t) is a 1+1 
dimensional field Q(x,t) subjected to the following chirality condition [43] 

^^- = d x Q(x l t). (34) 

A specific feature of the Hamiltonian formalism for non-local theories is that is 
contains the Euler-Lagrange equations as Hamiltonian constraints. 
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By using (34) we have 

(j t ) n Q(x,t) = (d x ) n Q(x,t),ne N , (35) 

and having in mind that Q(x, t) = q(x + t) assures the one-to-one correspondence 
between q(t) and Q(x,t) [43]. Ostrogradski's coordinates are defined as follows: 

Q in) (t) = (d x ) n Q(x,t)\ x=X0 , (36) 

where the discontinuity curve x (t) = x is a constant [43]. 

By using the inverse relation provided by the Taylor expansion around x = x 
we obtain 

QM = f j ^—^Q^\t). (37) 

n=0 Ul 

We notice that we reobtain the chirality condition as 

Q<ri(t) = Q ( - n+1 \t). (38) 

A boundary Poisson bracket was introduced in [43] and it has the expression 
as follows 



{F(t),G(t)}= E c k ,i iy4d x ) k+l [^§(x,t) S -^(x,t)]-(F^G), (39) 



k,l=0 



6QW 



where the coefficients Ck,i are constants and normalized in such a manner to satisfy 
the Jacobi identity. Fixing the time t the classical canonical relation is naturally 
obtained from (39) as 

{Q(x,t),P(x,t)} = 6 R (x-x). (40) 
By using (39) and (40) and defining Ostrogradski's momenta P( n )(t) as 

p (B) (t)= n dx {x ~ x ° )n p{x,t). (4i) 

J — oo Til 

From (41) the form of P(x, t) is given by 

oo 

P(x,t) = E P { n)(t)(-d x ) n 5 R (x - x ). (42) 

n=0 

By using (42) the expression for P( n ) (t) is as follows 

^(^E(- a .r-««|. (43, 
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The above definitions leads us to the following Hamilton equations 

P(n)(t) + P (B _D(t) = AT, (44) 

. g4[g]g) 

Taking into account the definition of P(o)(t) we observe that (45) is the Euler- 
Lagrange equation corresponding to the lagrangian Lj [Q] (t) written in the com- 
pact form. Finally, the Hamiltonian corresponding to L'j is given below 

oo 

H=Y, P (n)(t)Q {n+1) (t)-L' f [Q](t). (46) 

n=0 

Let us consider the following particular case 

L' f = ±( a T>?q(t)) 2 + a T>?F( q (t)). (47) 

By using (9) we obtain 

By making in (47) the corresponding substitutions of q(t) into Q(x, i) we obtain 

4<«M ^ H^ W,l) ^l(*)fl^^ 0W(i ' l) ) 2 + ^ F(W ' l)) 



By using (49) the corresponding fractional Euler-Lagrange equation has the 
form 



(49) 



+ 



(t-q)- Q / (t-q)- fl ^ _^ , ^ f a ^ (t - a) 
r(l - a) \r(i - a 
(t- a)- Q dF(Q(t,x)) 
r(l-a) dQ{t,x) 

' / a\,, {t-a) k ~ a k dF^{Q{t,x))^JL 1 /QW(t,a;) V* 



+ £U J* ! r(*-a + l)S dQ(f,J) E S 



r=1 Or! \ r! 



(50) 
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In the compact form (50) becomes 

(t-a)- dF(Q(t,x)) 
P(o)(t) = W ^ a B t Q(t,x)+ a B t dQ{tx) . (51) 

The fractional Hamiltonian equations are given below 



/a \ (* - a)' 

X 1 



\n J r(n - a + 1) 
°° / «\ ,, (t - a) k ~ a k 



SI k J fc! r(ife-a + l)§ rfQW(t,x) ^ 



fc=i 

and 



gW(t) = g(«+D(t). (53) 

We observe that (52) can be written in the compact form as 



P {n) (t) +P {n - 1) (t) = ( a r>?Q(t,x))i " • {l " y ' 




T(n-a + l) 



dQW(t,x) v 7 

We stress on the fact that x is not related with space [43]. The negative and the 
positive values of x corresponds to interactions with the past and respectively the 
future [43]. The solution of the chirality condition (53), will give Q(t, x) = q(x+t), 
which shows the explicit correspondence between q(t) and the left-mover Q(x, t). 

Taking into account that, when a — 1, T(l) = l,and using the fact that 

^ ^ ^ =0, for n > l,aDf(Q(t,x)) = d ^' t - > and the one-to one correspon- 
dence between q(t) and the left mover Q(x,t) we conclude that when the term 
aD^F(Q(t, x)) is switched off from (51) and (54), the classical results are reob- 
tained. 



5 Conclusions 

The fractional variational principles are powerful tools used successfully in various 
area of science and engineering. During the last years several points of view were 
launched in order to find the fractional Euler-Lagrange equations and to find an 
appropriate fractional Hamiltonian. One of the key point of these formulations 
was that under a certain limit the fractional theory includes the classical. As it 
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is expected the dynamics of the fractional calculus systems is different from the 
classical one but the classical dynamics is recovered as a particular case. 

In this paper we investigated the Euler-Lagrange and the fractional Hamil- 
ton equations corresponding to a fractional generalization of the equivalent La- 
grangians. Namely, we add to a given fractional Lagrangian a term which under 
certain limit reproduces the total derivative at the classical level. We have ob- 
served that for the fractional discrete systems the corresponding generalization of 
the classical equivalent Lagrangians lead us to a theory possessing infinite higher 
order derivatives. We have calculated the fractional Euler-Lagrangian for a par- 
ticular case as well as the fractional Hamilton equations were obtained in a most 
general case. 

We consider two specific examples in this paper. For the first example we 
consider one of the possible fractional generalization of the free one dimensional 
particle Lagrangian and we add to it a term of the form C a D"g(i). From the 
classical point of view the above investigated Lagrangians are equivalent. Both 
fractional Euler-Lagrange and Hamilton equations are obtained for this example 
and it was proved that they are equivalent in the fractional case. For example, 
the difference between the classical case and the fractional one is illustrated by 
the form of the fractional Euler-Lagrange equations given by (27). The second 
example deals with the same fractional generalization of the free one dimensional 
particle but we add a term of a D"F(g(t)). In this case, due to the Faa di Bruno 
formula, we used the 1+1 field formalism in order to obtain the corresponding 
fractional Euler-Lagrange and Hamilton equations. The classical results are re- 
obtained when a — 1. 
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